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In this article, we find out some analytical and numerical solutions to the problem of 
barrier tunneling for cluster deuterium, in particular using Langevin method to solve the 
time-independent Schrodinger equation. 



1 Introduction 

One of the most reported problem related to the CMNS (con- 
densed matter nuclear science, or LENR), is the low probabil- 
ity of Coulomb barrier tunneling. It is supposed by standard 
physics that tunneling is only possible at high enough energy 
(by solving Gamow function). 

However, a recent study by Takahashi (2008, 2009) and 
experiment by Arata etc. (2008) seem to suggest that it is 
not impossible to achieve a working experiment to create the 
CMNS process. 

In accordance with Takahashi’s EQPET/TSC model 
[1-3], the proposed study will find out some analytical and 
numerical solutions to the problem of barrier tunneling for 
cluster deuterium, in particular using Langevin method to 
solve the time-independent Schrodinger equation. It is hoped 
that the result can answer some of these mysteries. 

One of the results of recent experiments is the lack of 
signature of D-D reaction as in standard fusion process; this 
is part of the reason to suggest that D-D fusion doesn’t take 
place [1]. However, Takahashi suggests new possible reaction 
in the context of cluster deuterium, called 4D fusion [1-3], 
this mechanism seems to enable reaction at low temperature 
(CMNS). His result (2009) can be summarized as follows: 

“The ultimate condensation is possible only when the 
double Platonic symmetry of 4D/TSC is kept in its dy- 
namic motion. The sufficient increase (super screen- 
ing) of barrier factor is also only possible as far as the 
Platonic symmetric 4D/TSC system is kept. Therefore, 
there should be always 4 deuterons in barrier penetra- 
tion and fusion process, so that 4d simultaneous fusion 
should take place predominantly. The portion of 2D 
(usual) fusion rate is considered to be negligible”. 

In this respect it can be noted that there are recent reports 
suggesting that hydrogen cluster can get reaction at very low 
temperature, forming the condition of superfluidity [4]. This 
seems to happen too in the context of Takahashi TSC conden- 
sate dynamics. Other study worth mentioning here is one that 
discussed molecular chessboard dynamics of deuterium [5]. 

The difference between this proposed study and recent 
work of Takahashi based on Langevin equation for cluster 
deuterium is that we focus on solution of Schrodinger- 



Langevin equation [6, 7] with PT-Symmetric periodic poten- 
tial as we discussed in the preceding paper and its Gamow 
integral. The particular implications of this study to deuteron 
cluster will be discussed later. 

Another differing part from the previous study is that in 
this study we will also seek clues on possibility to consider 
this low probability problem as an example of self-organized 
criticality phenomena. In other words, the time required be- 
fore CMNS process can be observed is actually the time re- 
quired to trigger the critical phenomena. To our present 
knowledge, this kind of approach has never been studied be- 
fore, although self-organized criticality related to Schrodin- 
ger equation approximation to Burger’s turbulence has been 
discussed in Boldyrev [8]. Nonetheless there is recent study 
suggesting link between diffusion process and the self-orga- 
nized criticality phenomena. 

The result of this study will be useful to better understand- 
ing of anomalous phenomena behind Condensed matter nu- 
clear science. 

2 Schrodinger-Langevin equation 

The Langevin equation is considered as equivalent and there- 
fore has often been used to solve the time-independent Schro- 
dinger, in particular to study molecular dynamics. 

Here we only cite the known Langevin equation [3, p. 29] 

dX, — p,dt , (1) 

dp = -d x A 0 (X,)dt + Kp,dt + dW, y[7JK . (2) 

Takahashi and Yabuuchi also used quite similar form of 
the stochastic non-linear Langevin equation [7] in order to 
study the dynamics of TSC condensate motion. 

3 Schrodinger equation with PT-synunetric periodic po- 
tential 

Consider a PT-Symmetric potential of the form [9, 10] 

V = K\ sin(br ) , (3) 

where 
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Hence, the respective Schrodinger equation with this po- 
tential can be written as follows 

x ¥"(r) = — P(r) 'F(r) , (5) 

where 

2m 2m r n 

k(r) = —[E - V(r)] = -\E-k\ sm(fer)] . (6) 



The interpretation of this Gamow function is the tunnel- 
ing rate of the fusion reaction of cluster of deuterium (for the 
given data) corresponding to Takahashi data [12], with the 
difference that here we consider a PT-symmetric periodic po- 
tential. 

The numerical study will be performed with standard 
package like Maxima etc. Some plausible implications in 
cosmology modeling should also be discussed in the future. 



For the purpose of finding Gamow function, in area near 
x=a we can choose linear approximation for Coulomb poten- 
tial, such that 

V(x) - E = - a(x - a) . (7) 

Substitution to Schrodinger equation yields 

V" + 2 ff(x-a)V = 0, ( 8 ) 

Tv L 

which can be solved by virtue of Airy function. 



4 Gamow integral 

In principle, the Gamow function can be derived as foll- 
ows [11] 

d~y 

-^+P(x)y = 0. (9) 

Separating the variables and integrating, yields 



f ^- = j-P(x)Jx 



(10) 



or 



ydy — exp 




P(x) dx + C 



(11) 



To find solution of Gamow function, therefore the integral 
below must be evaluated: 



r = 






( 12 ) 
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For the purpose of analysis we use the same data from 
Takahashi’s EQPET model [3, 12], i.e. b = 5.6 fm, and r o = 5 
fm. Here we assume that E = V/, = 0.257 MeV. Therefore 
the integral becomes 



T = 0.218 Vm fk\ sin(for) - 0.257 dr . (13) 

Jr„ 

By setting boundary condition (either one or more of 
these conditions) 

(a) at r — 0 then Vo = — V* - 0.257MeV; 

(b) at r - 5.6 fm then V\ = k\ sin(/?r)- 0.257 = 0.257MeV, 
therefore, one can find estimate of m\ 

(c) Using this procedure solution of the equation (11) can 
be found. 
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